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> ■ ABSTRACT 

OO 



We describe a method for determination of physical conditions in the broad 
hne regions of a significant subsample of Seyfert-1 nuclei and quasars. Sev- 

O ' 6ral diagnostic ratios based on intermediate (AI111AI86O, Silll]A1892) and high 

(ClvA1549, SilvA1397) ionization hnes in the UV spectra of quasars are used to 
constrain density, ionization and metallicity of the emitting gas. We apply the 

L^ method to two extreme Population A quasars - the prototypical NLSyl I Zw 

5$ ; 1 and a high-z NLSyl-like object, SDSS J120144.36+011611.6. We find well- 

defined physical conditions: low ionization (ionization parameter < 10^^), high 
density (10^^ — 10^^ cm~'^) and significant metal enrichment. Ionization parame- 
ter and density can be derived independently for each source with an uncertainty 
that is always less than ±0.3 in logarithm. We use the product density times 
ionization parameter to estimate the broad line region radius and the virial black 
hole mass. Estimates of black hole masses based on the "photoionization" anal- 
ysis described in this paper are probably more accurate than those derived from 
the radius - luminosity correlation. 

Subject headings: quasars: general — quasars: individual (I Zw 1, SDSS J120144. 36+011611. 6) 
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Introduction 



We lack a simple well-defined diagnostic measure of physical conditions (density, ioniza- 
tion parameter, metallicity) in the broad line region of quasars. Techniques for est imating 



electron density and ionization in nebular astrophysics ( lOsterbrock fc Ferlandll2006[ ) are not 
straightforwardly applicable to the broad lines of quasars. Reasons include: hues are broad, 
suitable candidates are too closely spaced in wavelength, and density is at least an order 
of magnitude higher than the critical density for the forbidden transitions used in spectra 
of planetary nebulae and Hll regions. The ionization parameter can be estimated from the 
intensity ratio of two strong resonance lines arising from different ionic stages of the same 
element. It is again not easy to interpret the results in the case of quasars; for example, 
where the ratio Clll]A1909/ClvA1549 may not yield a meaningful value if the relatively low 
Clll]A1909 critical density implies that the the line is formed at larger radii than ClvA1549. 
Using lines of different ions introduces additional sources of uncertainty but widens the choice 
of diagnostic line ratios. 

Strong Clll]A1909 emission implies that He cannot be very high (ne~ 10^^~^'^cm~^). 
However, very high density is invoked to explain the rich low ionization spectrum (especially 
Fell) seen in quasars. Several lines in the UV spectrum of I Zw 1 (prominent Fell, relatively 
stron g AI111AI86O, and C l lllll76) point to wards high density at least for the LIL emitting 



zone (JBaldwin et al.lll996t iLaor et al.lll997l ). This low ionization BL R (or LIL BLR) has ver y 



similar properties to the Ol and Call emitting region identified by iMatsuoka et al.l ( 120081 ). 
The region where these low- ionization lines are produced cannot emit much Clll]A1909. But 
is Clll]A1909 really so strong in most quasars? BLR conditions are certainly complex and a 
single emitting region is not sufficient to explain both LILs and high ionization lines (HILs). 

In this letter we report on an analysis based on several diagnostic ratios used to con- 
strain density, ionization parameter, and metallicity in the BLR of two sources that are 
representative of relatively large sub-samples of quasars (©. These sources have the advan- 
tage to show weak or no Clll]A1909 greatly simplifying the interpretation of the emission 
line spectrum {^. Diagnostic ratios are heuristically defined (§!]) and interpreted through 
an array of photo-ionization simulations. We show that they converge toward well-defined 
values of ionization and density (§5]). The present analysis can be used to determine the 
product density times ionization parameter and the distance of the emitting region from the 
central continuum source tblr and the virial mass {^. 
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The sources 



The NLSyl I Zw 1 is an exam ple of the A3 spectral type (~ 10% of all quasars in the 
low-z sample of IZamfir et al.ll2010l ). Sources like I Zw 1 are also found at intermediate to 
high redshift. SDSS J120144. 36+011611. 6 at z ^3 is likely a high-redshift, high-luminosity 
replica of I Zw 1 with broader lines (Fig. [T]). Sources of spectral types A3 and A4 are 
the most extreme P opulation A quasars (FWHM of H/3 broad component ^ 4000 km s~^; 
Sulentic et al.l 120001 ) in terms of AI111AI86O and Fell emission. If we perform a search in the 
SDSS DR7 for quasars in the redshift range where both Civ A1549 and the A1900 blend 



are observed at optical wavelengths (2 



3) we find more than 200 sources out of 3000 



candidates with spectra resembling that of I Zw 1 on the basis of AI111AI86O strength. 



3. Analysis of the Emission Line Spectrum 

Felll lines are common and strong in the vicinit y of Clll]A1909 as seen in the SDSS 
template quasar spectrur n (IVanden Berk et al.l 120011) . They appear to be strong when 
AI111AI86O is also strong (jHartig fc Bald wini 1 19861 ). We adopt the template (option B) of 
Vestergaard fc Wilked ( 200ll ) fo r modeling Fell I emiss ion in our sources. Lya pumping 
enhances Felll (UV 34)11914.0 ([Johansson et al.l I2OOOI ) and this line can be a major con- 
tributor to the blend on the red side of Clll]A1909. The spectrum of I Zw 1 convincingly 
demonstrates this effect: both Clll]A1909 abd Felll (UV 34)11914.0 are needed to account 
for the double peaked feature at 1910 A that is too broad to be explained by a single line 
(Fig. [T])|^ The UV Fell template we adopt (Fell is not very strong in the spectral range we 
analyze) is based on a suitable CLOUDY simulation. 



The SPECFIT IRAF task (JKrissI 1 19941 ) allows us to fit the continuum, emission and 
absorption line components, as well as Fell and Felll contributions. We used the same 
power-law to describe the continuum at both the Civ Al 549 and 1900 A regions. We checked 
that the specfit solution properly accounted for the Fell UV 191 multiplet and the Felll 
features at 1930 and 2080 A that are unblended with other lines. 

The SPECFIT routine permits the deconvolution of Sim] A1892, making the AllllA1860/Silll] A1892 
ratio rehable. In Pop. A low (SillA1814) and intermediate ionization (AI111AI86O, Sim] A1892) 
lines (IILs) offer the key advantage of showing only the core broad component (BC) asso- 
ciated to low-ionization emission; HILs (ClvA1549 and SilvA1397) are decomposed into a 



^Note that the Galactic hne of Mg I at 2026 A can contribute to the spht appearance of the redshifted 
A1910 A feature (Fig. [1]). This absorption hne is included in our specfit analysis. 
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Lorentzian BC with the same width and shi fts of the IILs plus a b lueshifted residual (assumed 
to be Gaussian in the specfit procedure; iMarziani et al.ll2010l ). 



A second key advantage of the extreme Pop. A sources is that Cm] A1909 is weak or even 
absent (Fig. [1]). The fit to SDSS J120144.36+011611.6 suggests almost no Clll]A1909 emis- 
sion if enhanced Felll 11914 is included. If Felll 11914 is removed a maximum Cm] A1909 con- 
tribution is ^ 0.4 Silll]A1892 and the fit cannot model the red side of the blend. Clll]A1909 
emission in I Zw 1 is about ~ 0.6 Silll]A1892. 

The very dense region emitting the LILs should produce no Clll]A1909; so any emission 
from this line should arise in a different region. Reverberation mapping data for NGC 
5548 sugg est that Clll]A1909 i s mainly emitted farther from the central continuum source 
than H/3 (jPeterson et al.ll2004j ). There are also plausible physical scenarios th at can give 
rise to Clll]A1909 in a shielded, non- uniform environment (JMaiolino et al.ll2010l ). Since we 
are considering profiles integrated over the whole unresolved emitting region, the relevant 
question is: how much will the Clll]A1909 emitting regions contribute to the lines we are 
considering? The strongest Clll]A1909 emi tters along the El sequence (spectral type Al) 
typically show Silll]A1892/Clll]A1909 ^ 0.4 lKuraszkiewicz et all J2004h . Other IILs will be 
less affected. If the ionization parameter is logU ^ — 2.5, there will be also modest HIL 
emission. What we have estimated here is in many ways the maximum Clll]A1909 emission 
that we find weak. We therefore neglect the effect of any Clll]A1909 emitting gas on the 
IILs. 



4. Diagnostic ratios and their interpretation 



Figure [2] (left panel) shows the ionic fraction as a function of the geometrical depth in 
a slab of fixed column density (A^c= 10^^ cm~^) and density (nH= lO^^'^cm"^) exposed to a 
standard quasar continuum. A1+"'', Fe^"*", and Si^"^ share a region of dominance deep in the 
cloud, close to the end of the Stromgren sphere where the HILs are also emitted. It is very 
appropriate to consider IIL ratios since IILs are emitted in the same zones within the slab 
(right panel of Figure |2]) . 

Line ratios like AllllA1860/Silll]A1892 are useful diagnostics over a range of density that 
depends on their transition probabilities. Emission lines originating from forbidden or semi- 
forbidden transitions become collisionally quenched above the critical density and, hence, 
weaker than lines for which collisional effects are still negligible. The AllllA1860/Silll]A1892 
ratio is well suited to sample the density range 10^^ — 10^'^ cm~^. This corresponds to the 
densest, low ionization emitting regions likely associated with the production of Fell. 
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The ratios Sill]A1814 /Silll]A1892 and SilvA1397/ Silll]A1892 are sensitive to ionization, 
but roughly independent of metallicity. Conversely the ratio Civ A1549/Siiv A1397 is mainly 
sensitive to [Si/C]. 

We computed a multidimensional grid of C LOUDY simulations to derive estimate s 
of U and riu from our spectral measurements (JFerland et al.l Il998l : iKorista et al.l 119971 ) . 
CLOUDY computes population levels of the relevant ionic stages of Si, C, Al and, es- 
pecially for hues emitted in the fully ionized part of the slab (Fig. [2]), is expected to predict 
reliable relative line intensities. Simulations span the density range 7.00 < log nH< 14.00, 
and —4.50 < logU < 00.00, in intervals of 0.25 assuming of plane-parallel geometry. We 
considered three chemical compositions: (1) solar metallicity: (2) constant solar abundance 
ratio Al:Si:C with Z = 5Zq (5Z); (3) an overabundance of Si and Al with respect to carbon 
by a factor 3, again wit h Z = 5Zp^ (5ZSA). This last condition comes from the yields listed 
for type II Supernovae (IWoosley fc Weaven 119951 ). 



Results 



Tabled] reports the fluxes of the BC for IILs and HILs. We compute the ratios analyzed 
above from these values. They are shown as constant value contours in the U vs. nn plane 
obtained from cloudy simulations. At first, we assume metallicity equal to solar (Fig. El 
left panels). There is convergence towards a low ionization plus high density range. 

The discrepancy in the intersection point is significant for the SilvA1397/ClvA1549 ratio 
that depends mainly on the Si abundance relative to C. 

In equivalent plots made for Z = 5Z0SA, the agreement improves (Fig. [3]). High 
metallicity yields higher U and smaller nn if emission line ratios involving ClvA1549 are 
considered. This reflects the increase in abundance of Si and Al relative to C, and the fact 
the SillA1814, Silll]A1892, AI111AI86O lines are emitted at lower ionization than ClvA1549. 
The Si overabundance is also expected in the chemical composition of gas returned to the in- 
terstellar medium by a n evolved population w ith a top-loaded initial mass function simulated 
using Starburst 99 (JLeitherer et al.lll999l ). Independent evidence based on the strength 
of NvA1240 relative to ClvA1549 and Hel lA1640 indicates that chemical abundances may 
be 5 to 10 times solar (iDhanda et al.l 120071) in high redshift quasars, with Z ^ 5Zq reputed 
typical of high z quasars (JFerland et al.lll996l ). 



The halftone bands also show that we underestimate SillA1814 somewhat. SillA1814 
is the only low-ionization line considered in this study, and its formation is sensitive to the 
assumed X-ray continuum. IILs and HILs bypass the problem of LIL formation and offer 
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results consistent with a well-defined value of U and nn- 



6. Implications 

The ionization parameter U and tt-h are related to tblr through the equation tblr = 
const (52(t/?iH)^2 where Q is the number of H-ionizing photons. The dependence of U on 
^BLR was used to derive central black hole masses assuming a plausible a verage value of 



the product ni^-U and us ing the H/3 width as a virial broadening estimator (jPadovani et al. 



19901 : IWandel et al.lll999h . 



Using our derivations for nn and U, we obtain tblr and Mbh values li s ted in Table 
[21 Asymmetric errors have been quadratically propagated following iBarlowl (120041 ) . The 
last columns give the virial black hole masses follow ing our method and using the tblr ^ 
luminosity correlation ( IVestergaard fc Peterson! |2006| ) . The large intrinsic dispersion of the 
tblr ^ luminosity relation (±0.66 at a 2a confidence level) and its uncertain extrapolation 
makes the "photoionization" method presented here more appropriate for estimating Mbh 
in high-z quasars. 



7. CONCLUSION 

Diagnostics for density ionization and metallicity is possible for NLSyl-like sources at 
high redshift. Accurate diagnostics require high S/N and moderate dispersion but in principle 
can be applied to very high z (> 6.5) using data from IR spectrometers. The product {nnU) 
yields the possibility of deriving tblr and Mbh for a significant number of sources. Negrete et 
al. (2011, in preparation) will present an analysis of the applicability of the photoionization 
method described here to the general population of quasars. 
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Table 1. Measured quantities 



Object z /a(1700 A)^SiivA1397''CivA1549''SinA1814'"AliiiA1860''Siiii]A1892''Ciii]A1909'' FWHM<= 

I Zw 1 0.0605 3.1 19.3±4.0 27.0J;3;° S.st^;^ 16.6±3.3 31.1±6.2 18.0±3.6 1050±100 

SDSS J1201+0116 3.2332 1.6 12.9±2.6 15. Ol*'^ 3.0l^° 11.2±2.2 11.0±2.2 3.0: 4000±400 

''Rest-frame specific fiux in units of 10"^'* erg s~^ cm~^ A~^. 

Rest-frame flux in units of 10^^^ erg s~^ cm^'^. 
=Rest frame FWHM of IIL line BC and of the ClvA1549 BC in units of km s"!. 



Table 2. Derived quantities 



Object lognn'' log (7 log(nHl/) logrsLa'' logMfiH'^ logA/BH(VP06)d 

I Zw 1 12.05tg;Jg -2.65lg;^^ 9-4to;3 VI .2%^%^ 7-09to;2o 6.70 

SDSS J1201+0116 12.75lg;55 -3-00lg;ig 9.75lE!:^ 18.37l!5;i^ S-SStoi? 9-29 

''nn in units of cm~'^. 

''^BLR in units of cm. 

'^AfBH in units of M© computed using the customary virial formula Mbh= /^blrFWHM^/G, with the FWHM 
values of Tab. [T] and / = 0.75. 

'^A/bh in units of Mq computed following I Vestergaard fc PetersonI 1120061) . and input parameters reported in note 
(b). 
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Fig. 1.— Civ A 1549 (left) and A1900 spectral region (right) of I Zw 1 (top) and SDSS 
J120144. 36+011611. 6. The vertical scale is rest-frame specific flux in units of ergs s~^ cm~^ 
A~^ The major constituents of the blend are AhllA1860, Silll]A1892, Felll. The thin olive 
line represents template Fell emission, the thick one Felll. The dotted line corresponds to 
the Clll]A1909 estimated contribution. See text for details. 
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Fig. 2. — Left: Ionic Fractions as a function of the geometric depth h in the gas slab. Right: 
Line emissivity per unit volume in units of ergs s~^ cm~^ multiplied by depth h in cm. 
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Fig. 3. — Iso-contour of diagnostic ratios for I Zw 1 (top) and SDSS J120144. 36+011611. 6 
(bottom). The left panels refer to the case of solar metallicity; the right ones to the case 
of five times solar plus a factor of three metal enrichment (5ZSA). The grey band is the 
uncertainty band of the ratio SillA1814/Silll]A1892. Note that the SilvA1397/ClvA1549 is 
not proving a useful constraint in the 5ZSA case since it varies little in a wide area around 
the intersection point in the (f/,nH) plane. 



